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Abstract 

We  present  a  theoretical  model  of  the  coupled  evolution  of  the  phase  space  densities  of  solar  energetic 
particles  (SEP)  and  interplanetary  Alfven  wave  distribution  for  interpreting  the  evolution  of  elemental 
abundances  in  SEP  accelerated  by  a  traveling  interplanetary  shock.  The  model  is  formulated  in 
“mixed”  coordinates  and  includes  the  effects  of  magnetic  focusing,  pitch-angle  scattering,  solar-wind 
convection,  and  adiabatic  deceleration,  as  well  as  wave  transport  and  amplification  by  streaming 
protons.  The  results  suggest  that  a  key  ingredient  for  understanding  SEP  abundance  variation  is  the 
transport  of  shock  accelerated  ions  through  proton  excited  non-Kolmogorov  wave  spectra. 


1  Introduction 

The  acceleration  and  transport  of  SEP  has  long  been  a  challenge  for  theoretical  modeling  because 
of  the  variety  and  complexity  of  the  physical  processes  involved  (Lee,  1997;  Reames,  1999).  Despite 
the  diffculties,  the  challenge  must  be  met  not  only  for  theoretical  interest  but  also  for  interpreting  SEP 
observations.  The  pressing  need  for  a  realistic  model  is  manifest  when  we  attempt  to  understand  the 
systematic  variations  of  SEP  abundances  observed  by  the  sensitive  Wind/EPACT  instrument  (Tylka, 
Reames  &  Ng,  1999a, b).  The  temporal  variations  of  SEP  abundances  explore  the  differences  in  the 
acceleration  and  transport  experienced  by  SEP  ions  of  different  A/Q  ratios  (A  =  atomic  mass,  Q  = 
charge  state),  and  challenge  our  ingenuity  in  separating  the  roles  played  by  source  abundances,  charge 
states,  acceleration  and  interplanetary  transport.  We  have  taken  two  important  steps  in  modeling  the 
evolution  of  ion  phase  space  densities  by  including  (i)  the  effect  of  self  excited  waves  on  SEP  transport 
(Ng  &  Reames,  1994),  and  (ii)  a  moving  source  of  energetic  ions  associated  with  a  traveling  shock 
(Ng,  Reames  &  Tylka,  1999).  In  this  paper,  we  present  preliminary  results  from  our  newest  model 
which  now  includes  the  effects  of  solar-wind  convection  on  ions  and  waves,  and  adiabatic  deceleration. 
The  new  results  support  our  earlier  suggestion  that  ion  transport  through  proton  amplified  waves  is 
key  to  understanding  the  temporal  evolution  of  SEP  abundances. 

2  Model  and  Equations 

The  model  includes  magnetic  focusing,  pitch-angle  scattering,  wave  amplification,  energetic  ion 
source  at  a  traveling  “shock” ,  solar-wind  convection  of  waves  and  ions,  and  adiabatic  deceleration.  We 
make  several  simplifications  to  render  the  model  tractable.  Acceleration  is  decoupled  from  transport 
by  injecting  time  varying  power-law  energetic  ion  spectra  at  a  “shock”  that  travels  with  constant 
velocity  Vjh  on  a  radial  magnetic  field.  The  phase-space  densities  fffp,  P,  r,  t)  of  ion  species  s  and  the 
differential  magnetic  energy  density  La(k ,  r,  t)  of  er-mode  Alfven  wave  satisfy  the  following  differential 
equations  in  “mixed  coordinates”:  (r,  t)  in  fixed  inertial  frame  and  (p,  P)  in  the  solar-wind  frame. 
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In  the  above,  t  is  time,  r  heliocentric  distance,  v  ion  velocity,  P  rigidity,  /i  pitch-angle  cosine, 
D'li/iiP-  fJ  - E  t)  species-s  diffusion  coefficient  in  //-space.  B  =  i?o(ro/'r)2  magnetic  field,  k  wavenumber, 
7]  =  In (k/B),  G s  ion  source  term,  7 a  growth  rate  of  Ia,  c  light  speed,  e  elementary  charge,  E  proton 
total  energy,  Va  oc  r  1  Alfven  speed,  W  solar-wind  speed,  I4h  shock  speed,  rs h  =  roiSh  +  tVs\A  shock 
distance  from  Sun,  Va  =  W  +  QoVa,  go  =  ±1  for  outward  (inward)  waves,  Bq,  Uq,  tq  and  ro,sh  are 
constants,  and  /?(/,(//,  v.  P.  k.  V„)  is  the  wave-particle  resonance  function  (Ng  &  Reames,  1995).  In 
eq.  (1),  the  third,  fourth  and  fifth  terms  represent  the  effects  of  focusing,  adiabatic  deceleration,  and 
pitch-angle  diffusion  respectively.  See  Ruffolo(1995),  Barnes  (1979)  for  further  discussion  on  eqs.  (1) 
and  (2),  respectively.  Only  /h  appears  in  eq.  (4)  since  the  minor  ions  contribute  negligibly  to  wave 
evolution  (Lee,  1983).  Eq.  (1)  for  /h  and  eq.  (2)  for  Ia  are  coupled  via  eqs.  (3)  and  (4).  The  ion 
source  Gs  per  (cm2  MV3  h)  is  assumed  isotropic  and  sharply  peaked  at  r  =  r^,  such  that 
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and  Gs  is  negligible  outside  (rsh  —  Ar/2,  rsh  +  Ar/2),  with  A r  the  finite-difference  grid  size.  In  eq.  (5), 
ns(ro)r(\/ >'2  is  the  solar-wind  s-ion  density  at  r,  r 0  =  1  AU,  PqjS  is  the  “seed”  ion  rigidity,  and  we 
set  a  =  8  x  10~5  h_1  for  all  ion  species.  We  take  uh(^o),  no(ro)  and  //pefco)  in  proportion  to  the 
observed  coronal  abundances  (Reames,  1995),  e.g.  6.4  cnr3,  4.1  x  ICE3  cnr3  and  5.5  x  ICE4  cnr3, 
respectively.  All  ions  are  seeded  at  the  same  5  keV/amu  (Lee,  1983). 

For  the  results  in  Section  3,  we  start  with  empty  ion  distributions  and  forward  propagating  right 
and  left-hand  polarized  Alfven  wave  intensities  which  are  steady-state  solutions  of  eq.(2)  with  Kol¬ 
mogorov  spectra,  such  that  the  ion  mean  free  paths  A  =  3/8  dg(l— g2)2v / >  1  AU  in  the  inner 
heliosphere.  At  t  =  0,  a  “shock”  travels  out  at  specified  velocity  V^h,  releasing  energetic  ions  with 
spectral  index  a  =  6  +  f3Vsht  (eq.[5]).  Classical  shock  acceleration  theory  gives  a  =  —3 cr/(cr  —  1), 
with  cr  the  compression  ratio,  and  a  =  —4  for  the  maximum  cr  =  4.  The  coupled  eqs.  (1)  and  (2) 
are  solved  using  a  refinement  of  the  finite-difference  techniques  described  in  Ng  &  Reames  (1994).  In 
Section  3  we  study  how  the  Fe/O  and  He/H  ratios  depend  on  the  shock-related  source  strength. 


3  Results 

We  study  the  evolution  of  Fe/O  at  2.6,  3.7,  and  5.2  MeV /amu,  and  He/H  at  2.2  MeV/amu,  in  the 
energy  range  reported  by  Tylka  et  al.  (1999a, b).  We  take  Ape/QFe  =  4.  Ao/Qo  =  Ane/Qi-ie  =  2. 
The  calculated  abundance  ratios  are  displayed  versus  time  in  Fig.  l(a)-(d).  In  all  cases,  all  Fe/O 
ratios  decrease  precipitously  after  event  onset.  For  moderately  strong  but  fast  weakening  shock  (a), 
they  continue  to  decrease  at  slower  rates,  punctuated  by  a  brief  small  “shock  spike”.  For  maximum 
compression  shock  that  weakens  rapidly  (b),  all  Fe/O  ratios  rebound  to  energy-dependent  values 
>  2,  then  decrease  to  <  1,  faster  at  higher  energy.  For  the  persistent  maximum  compression  shock 
(c),  all  Fe/O  ratios  rebound  to  sa  2,  decaying  very  slowly  past  unity  at  shock  passage.  The  energy 
dependence  of  Fe/O  is  different  in  case  (c)  than  cases  (a)  and  (b). 

At  onset,  He/H  decreases  rapidly  like  Fe/O  for  the  weak  shock  in  Fig.  1(a),  but  rises  sharply  for 
the  stronger  shocks  in  (b)  and  (c).  Since  at  the  same  MeV/amu,  He2+  (Fe14+)  has  higher  rigidity 
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Figure  1:  Fe/O  and  He/H  ratios  versus  time  at  r  =  1.125  AU  for  different  source  strength 
of  shock-accelerated  ions,  (a)  Initially  moderately  strong  but  fast  weakening  shock, 
(b)  Initially  maximum  compression  but  also  fast  weakening  shock,  (c)  Persistent 
maximum  compression  shock.  The  ion  source  momentum  power-law  index  is  given  by 
a-  =  8  +  /3Vsh  t  when  the  shock  is  at  rsh  =  roiSh  +  il 4h-  (d)  Similar  to  (b),  but  with  no 
solar- wind  convection  or  wave  transport.  Vertical  arrow  indicates  shock  passage. 
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Figure  2:  (a)  Differential  magnetic  intensity  of  forward  right-hand  Alfven  wave  at  r  = 

0.35 AU.  (b)  Evolution  of  radial  profiles  of  mean  free  path  of  2.6  MeV/amu  Os+. 

and  longer  mean  free  path  than  H+  (Os+)  in  a  Kolmogorov  wave  spectrum,  this  behavior  is  at  first 
sight  puzzling.  In  fact,  both  types  of  He/H  variation  at  onset  (in  contrast  to  the  decreasing  Fe/O) 
have  been  reported  (Tvlka  et  ah,  1999a, b;  Mason  et  al.,  1983,  Witte  et  ah,  1979). 

For  comparison,  we  display  in  Fig.  1(d)  the  results  obtained  without  solar-wind  convection  and 
wave  transport,  with  parameters  similar  to  those  in  Fig.  1(b).  The  two  sets  of  results  are  qualitatively 
similar.  However,  He/H  shows  in  Fig.  1(d)  a  valley  which  is  absent  in  Fig.  1(b).  Although  we  do  not 
display  them  here,  the  proton  intensity  time  profiles  predicted  by  the  model  in  these  calculations  are 
in  the  range  of  the  observations. 


Fig.  2(a)  shows  the  evolution  of  right-hand  polarized  Alfven  wave  spectrum  at  r  =  0.35  AU 
from  an  initial  Kolmogorov  spectrum,  while  Fig.  2(b)  shows  the  evolution  of  the  radial  profiles  of  2.6 
MeV/amu  oxygen  mean  free  path  to  t  =  36  h,  both  for  the  parameters  of  Fig.  1(b).  These  two  figures 
reveal  the  dynamic  response  of  Alfven  waves  to  energetic  proton  streaming,  and  argue  persuasively 
that  this  must  have  important  consequences  on  SEP  transport  and  abundance  evolution. 

4  Discussion 

The  key  to  understanding  the  three-stage  abundance  variations  (Fig.  1)  with  the  contrasting 
behavior  between  2.6  MeV/amu  Fe/O  and  2.2  MeV/amu  He/H  can  be  gleaned  from  Fig.  2,  viz. 
proton  excitation  of  Alfven  waves.  The  precipitous  fall  of  Fe/O  (and  He/H  if  it  does)  is  simply  due  to 
propagation  effect  in  the  background  Kolmogorov  wave  field  (oc  k.-5/3).  For  I a  oc  k  b,  the  quasi-linear 
theory  predicts  D^fv  oc  Pb~2.  Hence  D^/ji  increases  (decreases)  with  P  if  —  b  is  >  (<)  —  2  at  the 
resonant  wavenumber  kres.  A  k~2  spectrum  is  shown  in  Fig.  2(a)  for  comparison.  The  background 
fc~5/3  spectrum  implies  Dllfl/v  oc  P1/3  everywhere  and  so  Fe/O  and  He/H  will  fall  precipitously 
after  event  onset.  However,  Fig.  2(a)  shows  that  wave  evolution  brings  in  D^/ji  that  decreases 
with  P  at  some  kres  and  increases  with  P  at  other  kres,  in  a  temporally  (and  spatially)  dependent 
manner.  The  rise  of  He/H  may  thus  be  explained  as  follows.  Protons  with  twice  the  speed  and  the 
same  rigidity  as  2.2  MeV/amu  He2+  and  /i  ~  1  race  ahead  at  event  onset  to  amplify  an  mute  Alfven 
waves  that  subsequently  scatter  the  later  arriving  He2+.  With  sufficiently  high  proton  intensity  from 
a  strong  shock,  this  wave  growth  causes  the  first  “wave”  of  2.2  MeV/amu  He2+  to  arrive  later  than 
the  first  “wave”  of  2.2  MeV  protons,  and  so  He/H  rises  from  <  1.  However,  the  proton  driven  wave 
growth  is  not  large  enough  at  higher  rigidities  (lower  wavenumbers)  to  reverse  the  behavior  of  Fe/O 
soon  after  event  onset.  The  rise  of  Fe/O  in  the  second  phase  is  due  to  the  continuing  wave  growth 
retarding  Os+  more  than  Fe  l4+.  The  final  decline  as  the  shock  approaches  is  due  to  more  Os+  being 
retained  near/behind  the  shock  than  Fe14+.  The  later  He/H  variation  may  be  understood  similarly. 

We  will  continue  to  study  the  model  here  to  understand  its  full  implications  and  to  delineate  its 
range  of  validity.  We  will  also  extend  the  model  to  include  nonlinear  wave  cascading. 
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